Magnetic phase diagram of the frustrated S = 1/2 chain magnet LiCu2 02 
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We present the results of the magnetization and dielectric constant measurements on untwinned 
single crystal samples of the frustrated S = 1/2 chain cuprate LiCu2 02. Novel magnetic phase 
transitions were observed. A spin flop transition of the spiral spin plane was observed for the field 
orientations H \\ a, b. The second magnetic transition was observed at H « 15 T for all three 
principal field directions. This high field magnetic phase is discussed as a coUinear spin-modulated 
phase which is expected for an 5=1/2 nearest-neighbor ferromagnetic and next-nearest-neighbor 
antiferromagnetic chain system. 

PACS numbers: 75.50. Ee, 76.60.-k, 75.10. Jm, 75.10.Pq 



I. INTRODUCTION 

Unconventional magnetic orders and phases in frus- 
trated quantum-spin chains are one of the attractive is- 
sues, because they appear under a fine balance of the ex- 
change interactions and are sometimes caused by much 
weaker interactions or fluctuationsir— . 

A kind of frustration in quasi-one-dimensional (ID) 
chain magnets is provided by competing interactions 
when the nearest neighbor (NN) exchange is ferromag- 
netic and the next-nearest neighbor (NNN) exchange is 
antiferromagnetic. Numerical investigations of the frus- 
trated chain magnets within different models^"— have pre- 
dicted a number of exotic magnetic phases in the magne- 
tization process such as planar spiral and different multi- 
polar phases. Theoretical studies of the magnetic phase 
diagram show that the magnetic phases are very sensitive 
to the interchain interactions and anisotropic interactions 
in the systeroiSrii. 

There are a number of magnets which are at- 
tractive objects for the experimental investigations 
as realizations of one-dimensional frustrated sys- 
tems: LiCuV04, RbzCuzMoaOiz, Li2ZrCu04, CuClz, 
Li2Cu02, LiCu202, NaCu202 (see for example, RefslTB- 
[Toh . Some of these magnets, such as LiCuV04, Li2Cu02 
and Li2ZrCu04, are quasi-one-dimensional magnets, 
while others are quasi-two-dimensional magnets with the 
layers of frustrated chains. 

In the present paper, we report the results of magneti- 
zation and dielectric constant measurements on the quasi 
two dimensional antiferromagnet LiCu202. Our new ex- 
periments expand the magnetic field range from the fields 
of up to 9 T of the previous work^S to the fields of up to 
52 T, which is approximately 40% of the expected satura- 
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FIG. 1: (a) Schematic representation of the arrangement of 
Cu^^ moments in LiCu202. Only one of four Cu^^ (a, /3, 7, 
(5) positions is shown. Arrows correspond to the spin direc- 
tions below Ta2- J\,Ji and J3 are the main exchange inte- 
gralsi^. (b) Scheme of the orientation of the spin plane for 
H in the (6c)-plane. H > H^i- Spin plane is shadowed with 
gray. The direction of spin rotation by displacement along 
the chain (6-axes) is shown by the circular arrow, n is the 
vector normal to the spin plane. 



tion field {Haat ) ■ These measurements display two phase 
transitions. The low field transition may be interpreted 
as a spin-flop transition while the novel high field transi- 
tion is, most probably, of the exchange origin. Using the 
results of a ID frustrated model, we suggest that the ob- 
served high fleld transition from a spiral magnetic phase 
is the transition to a collinear spin-modulated magnetic 
phase. 
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FIG. 2: Upper panehPolarized optical microscope images of 
the (a6)-plane of LiCu2 02 samples (a) without twinning, (b) 
with two large twinned domains and (c) with complicated 
twinning structure. Lower panel: X-ray diffraction pattern 
observed in 6 — 26 geometry on untwinned single crystal. The 
spectra were obtained from the (be) and the (ac)-planes of 
the crystal. 



II. CRYSTALLOGRAPHIC AND MAGNETIC 
STRUCTURES OF LICU2O2 

LiCu202 crystallizes in the orthorhombic lattice (space 
group Pnma) with the unit cell parameters a =5.73 A, 
h =2.86 A and c =12.42 A.— The unit cell parameter a 
is approximately twice the unit cell parameter h. Con- 
sequently, the LiCu202 samples, as a rule, are charac- 
terized by the twinning due to formation of the crystal- 
lographic domains rotated by 90° around their common 
crystallographic axis c. 

The unit cell of the LiCu2 02 crystal contains four uni- 
valent nonmagnetic cations Cu"*" and four divalent cations 
Cu^+ with the spin 5*= 1/2. There are four equivalent 
crystallographic positions of the magnetic Cu^"'' ions in 
the crystal unit cell of LiCu2 02, conventionally denoted 
as a, 13, 7, and S. The chains formed by one of four kinds 
of Cu^^ ions and relevant exchange interactions^^ within 
the system are shown schematically in Fig. [IJa). 

The two-stage transition into a magnetically ordered 
state occurs at = 24.6 K and Tc2 = 23.2 Neu- 
tron scattering and NMR experiments have revealed an 
incommensurate magnetic structure in the magnetically 
ordered state (T < rci)^i2ii^. The wave vector of the 
incommensurate magnetic structure coincides with the 
chain direction (&-axis). The magnitude of the propa- 
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FIG. 3: (color online) Upper panel; Temperature dependence 
of the low field magnetization curves M{T) / fioH for three 
field directions H \\ a, b, c and noH = 0.1 T. Lower panel: 
M{T)/fioH in the vicinity of the transition temperatures Td 
and Tc2 at hqH =1, 3, 5, and 7 T for || b. 



gation vector at T < 17 K is almost temperature inde- 
pendent and is equal to 0.827x27r/6. The neutron scat- 
tering experiments have shown that the magnetic mo- 
ments neighboring along the a-direction are antiparal- 
Icl, whereas those neighboring along the c-direction are 
coaligned. The intra-chain and inter-chain exchange con- 
stants were determined from the analysis of the spin 
wave spectraJ^ The nearest neighbour exchange inter- 
action is ferromagnetic Ji = —7.00 mcV, while the next 
nearest neighbour exchange interaction is antiferromag- 
netic J2 = 3.75 meV. The competition between these in- 
teractions leads to incommensurate magnetic structure. 
The antiparallel orientation of the magnetic moments of 
Cu^+ of the neighbouring chains is caused by the strong 
antiferromagnetic exchange interaction J3 = 3.4 meV. 
The coupling of the Cu^"'' moments along the c-direction 
and the couplings between the magnetic ions in different 
crystallographic positions are much weaker— i^. Thus, 
LiCu202 can be considered as a quasi- two dimensional. 
The quasi-two-dimensional character of magnetic inter- 
actions in LiCu2 02 compound was also proved by the res- 
onant soft X-ray magnetic scattering experiments ^^1^^ . 
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FIG. 4: (color online) Temperature dependences of the mag- 
netization derivative dM/dT for three field directions in the 
vicinity of Tci and Tc2 obtained from M{T) data at ^.qH =1, 
5, and 7 T. 



Magnetic structure of LiCu202 at zero magnetic field 
was studied by several groups by means of neutron 
diffraction experiments22i22iiS. The authors of Ref. [2^ 
have proposed the planar spiral spin structure with the 
spins confined to the (a6)-plane. Polarized neutron scat- 
tering measurements of Reff23l have detected the spin 
component along the c~direction, indicating the spiral 



magnetic structure in the (6c)-plane. The authors of 
Ref. [25I alternatively, have proposed a spiral spin struc- 
ture confined to the (1,1,0) plane. It was attempted to 
extract information about magnetic structure of LiCu202 
from the studies of electric polarization, which accompa- 
nies the magnetic orderin g^^'^^ . Unfortunately, at the 
moment, the nature of this polarization is not clear—, 
which did not permit to obtain the definitive information 
about the zero-field magnetic structure from this type of 
experiment. 

The magnetic structure realized in LiCu2 02 is much 
clearer in the magnetic fields above 3 T applied along 
b and c directions. ESR and NMR studies of LiCu202 
in the low-temperature magnetically ordered phase (T < 
Tc2) have shown that the planar spiral magnetic structure 
is formed in this compoundSS. The magnetic moments 
located at the a (/3, 7 or 5) position of the crystal unit 
cell with coordinates x, y, z (measured along the a, b and 
c- axes of the crystal respectively) are defined as: 



Ma = /U • hi-l)'-^^" ■ COs(fcjc • y + 0a) + 

+//-Z2(-l)"/"-sin(fc,,. 2/ + </.„), (1) 

where li and I2 are the two mutually perpendicular unit 
vectors, kic is the incommensurability vector parallel to 
the chain direction (6-axis) and fi is the magnetic moment 
of the Cu^"*" ion. The ordered magnetic moment per cop- 
per ion at r < 10 K was evaluated as ^ = 0.85 /zb.— 1^ 



The phases 



and 6x determine the mutual ori- 



entation of the spins in the chains formed by the ions in 
the different crystallographic positions. Their values can 
be extracted from the NMR dataSSi^. 



III. SAMPLE PREPARATION AND 
EXPERIMENTAL DETAILS 

Single crystals of LiCu202 with the size of several cu- 
bic millimeters were prepared by the solution in the melt 
method as described in Ref. [l^. The samples had a shape 
of a flat plate with the developed (a6)-plane. The twin- 
ning structure of the samples was studied by the optical 
polarization microscopy. Typical views of the (a6)-plane 
of the samples are shown in Fig. [21 The dark and bright 
areas on the photos correspond to the domains with two 
different directions of the a-axis. In most cases, the sam- 
ples were twinned with characteristic domain size of sev- 
eral microns (see Fig.^^c)). It was possible to select the 
samples with large domains (see Fig. [2{b)) and the best 
samples were without twinning structure (Fig. [21[a)). 

The absence of twinning structure in the samples 
selected for the experiments was confirmed by X-ray 
diffraction measurement (see lower panel of Fig. [2]). The 
X-ray diffraction patterns were taken in 9 — 29 geometry 
with use of CuK(q;) irradiation. The untwinned crys- 
tal produced clearly distinguishable diffraction patterns 
from the (be) and the (ac)-planes of the crystal (see lower 
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FIG. 5: (color online) Upper panel: Magnetization curves 
M{H) at 4.2 K for three field directions H \\ a, b, c measured 
in pulse magnetic fields of up to 52 T. The inset shows the 
field dependence of dM/dH at 4.5 K for H \\ h measured with 
a SQUID magnetometer. Lower panels: Measured magnetic 
field dependences of dM/dH for the principal field directions 
H II a, b, c. The inset shows the enlarged dependency in the 
range of transition field Hc2 for H \\ b. 



panel of Fig. [2]). The spectra obtained from twinned sam- 
ples had a form of superposition of these spectra. ESR 
measurements also confirm absence of twinning in the 
samples used^S. 

Magnetization curves in static magnetic fields of up 
to 7 T were measured with a commercial SQUID mag- 
netometer (Quantum Design MPMS-XL7). High field 
magnetization data were taken in magnetic fields of up 
to 52 T using a pulsed magnet at KYOKUGEN in Osaka 
Universityi^ 

The dielectric constant e was obtained from the four- 
contact capacitance measurements at 10 kHz on the ca- 
pacitor formed by contact plates attached to the sample. 
The naturally grown samples with the fiat (a&)-planes 
were used for the measurements of Ec- Samples for the 
measurements of and e^, were properly cut from the 
single crystals without twinning structure. Capacitor 
plates were formed right on the sample surface either by 
the sputtering of the silver film, or with the silver paste. 

Different samples from the same batch were used for 
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FIG. 6: (color online) Upper panel: Temperature depen- 
dences of dynamic dielectric constants of three samples at 
zero magnetic field for alternating electric field E \\ a,c. De- 
pendences e{T) presented in the main panel were obtained 
by subtraction of the linear temperature drifts. Curves are 
shifted along the y-direction for better presentation. Insert 
shows the raw temperature dependence of the capacitance 
C(T) for £ II a at H = with the linear extrapolation of 
C(T) at T > Tc to the range of T < Tc2 shown by dotted 
line. Lower panel: The field dependences of the capacitance 
for II a and the magnetic field H \\ b a.t different temper- 
atures. The data in the lower panel are plotted without shift 
along the y-axis. 



the dielectric constant and magnetization measurements, 
these samples are further denoted as " sample 1-5" . 



IV. EXPERIMENTAL RESULTS 

A. Magnetization measurements at low fields 

The experimental curves M[T)/^,qH and their tem- 
perature derivatives for three field directions H \\ a,b, c 
at i^lqH = 0.1 T are shown in Figs. [3] and |31 For 
all three field directions, there is a broad maximum of 
M{T)/^qH at T = 38 K which is a typical feature of 
low-dimensional antifcrromagncts. In the paramagnetic 
region, magnetic susceptibility for H || c exceeds that 
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for H II a, b, which is consistent with the anisotropy of 
the p-tensor measured in ESR experiments^^ {ga,b = 2.0; 
Pc = 2.2). The magnetic susceptibihties are tempera- 
ture independent for T <C Tjv . The absence of significant 
paramagnetic contribution at low temperatures indicates 
low concentration of the uncontrolled magnetic impuri- 
ties. 

The transitions to the magnetically ordered states at 
Tci and Tc2 are marked by the change of the M{T) slope. 
These inflection points are well resolved on the temper- 
ature derivative of M/hqH (Fig. SI). The -^{M/^iqH) 
curves are strongly anisotropic. For if || b and c, there 
are two anomalies corresponding to two transitions at Td 
and Tc2- The temperature dependence of -^{M/ fioH) 
for II a shows only one sharp peak at the tempera- 
ture near Tc2 in all studied field range. An additional 
anomaly was observed on the M{T) curves for || b 
and /io-ff 3 T (lower panel of Fig. At the transition 
temperature Tc3, the slope of A/(T) changes abruptly. 
The value of Tcs increases with the increase of magnetic 
field and merges with Td at > 6 T. 
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FIG. 7: (color online) Magnetic phase diagrams of LiCu202 
compound for H \\ h (upper panel) and .ff |j c (lower 
panel). The different symbols in the figures correspond to 
the data from different experimental methods. Solid triangles 
and circles show the phase boundaries obtained from ec{T) 
and ea{H) measurements, respectively. Open red circles and 
squares show the phase boundaries obtained from M{H) and 
M{T) measurements, respectively. The solid lines are guides 
to the eye. The grey shadowed area in the upper panel marks 
sample-dependent position of the phase boundary between 
phases I and II. High field spin structure III schematically 
shown in the phase diagram is merely speculative. Dashed 
line is a schematic phase boundary between the paramagnetic 
polarized and the ordered phases. 



B. High- field magnetization measurements 

The upper panel of Fig. [5] demonstrates magnetization 
M{H) curves measured in pulsed fields of up to 52 T for 
H II a,b,c. These curves were obtained by integrating 
the dM/dH curves presented in the lower three panels of 
Fig. El 

The dM/dH curves show two anomalies which are 
marked as Hd and Hc2- The anomaly at Hd was ob- 
served for H \\ a and for H \\ h. It is characterized 
by the increase of the dM{H)/dH, which may indicate 
a spin-reorientation transition. At the second transition 
field Hc2, observed in all three principal orientations, the 
field derivative of magnetization decreases. 

The position of the critical field Hd for the J? || b cor- 
responds to the magnetic transition at the temperature 
Tc3 observed in the temperature scan of magnetization. 
The transition field Hd is strongly anisotropic: for the 
H II a, i-ioHd = 15 T, while for the i? || ft it varies from 
2 to 5 T from sample to sample and demonstrates a hys- 
teresis. The differential susceptibility dM/dH increases 
stepwise at the transition field Hd by 65±5 % for || a 
and by 10 ± 3% for H \\ b. 

The second transition field Hc2 is less anisotropic, its 
values are ijlqHc2 = 15.3,16.7 and 11.5 T for the H \\ 
a,b, c, respectively. At this transition the differential 
susceptibility decreases stepwise by 25 ±3 %, 3±1% and 
11 ± 2% for the H \\ a, b, c, respectively. 



C. Dielectric constant measurements 

Magnetic and electric properties of LiCu2 02 are 
strongly couplc d^^'^^ . The magnetic transitions are ac- 
companied by anomalies of the dielectric constant. 



6 



The temperature and field dependences of the dielec- 
tric constant measured in our experiments are shown in 
Fig.El These curves were obtained from the temperature 
dependences of the capacitance. Linear drift obtained by 
extrapolation of C [T] aiT > Tc was subtracted for bet- 
ter presentation. The example of the raw C (T) curve and 
its high temperature linear extrapolation is given on the 
insert to Fig. [51 The temperature drift depends on the 
quality of capacitor plates contact with the sample. The 
strongest temperature drifts were observed for samples 
prepared for measurements at £^ || a,b. In these cases 
the samples were composed from 10 or more oriented 
parts fastened with glue. 

There are clearly visible anomalies on e{T) a.t E \\ a, c 
while there is no anomaly at E \\ b. The edT) demon- 
strates two peaks at Td and Tc2 as shown in the upper 
panel of Fig. |6l These peaks correspond to the two suc- 
cessive transitions between the paramagnetic phase and 
the ordered spiral one. The anomaly at Td has no hys- 
teresis, while the edT) shows hysteresis near Tc2- The 
positions of the peaks of Sc{T) varies from the sample to 
the sample in the range of 24.7 ± 0.2 K and 22.7 ± 0.2 K 
for Tci and Tc2, respectively. The peak value of the di- 
electric constant changes by the factor of 10 depending 
on the sample. The £a{T) shows only one broad, step- 
like anomaly around Tc2- The positions of the peak of 
ec{T) and of the maximal slope of £a{T) coincide for the 
same sample. No anomalies on the temperature and field 
dependences of Eh were observed. 

One broad peak is observed in the field dependence of 
capacitance C of LiCu2 02 sample for H \\ b and || a 
(lower panel of Fig. [6l), which corresponds to the tran- 
sition at Tc3 and Hd on the corresponding M(T) and 
M[H) dependences (see Figs. [3] and [5]). 



D. Phase diagram 

The H — T magnetic phase diagram of LiCu2 02 is pre- 
sented in Fig. [T] In total, it contains 5 different phases: 
(I) the low-field phase, (II) the phase above Hd for 
H \\ b, (III) the phase above H^2 for H || c, (IV) the 
intermediate ordered phase between Td and Tc2, (V) the 
polarized paramagnetic phase. 

The anomalies in M{H), M{T), ea{T), ea{H) which 
indicate the transition from the phase I to the phase II on 
the phase diagram have a hysteretic character. The phase 
boundaries obtained for three different samples for i? || b 
are also shown on the phase diagram. The boundary 
between phases I and II for iJ || b changes with the 
sample. 



V. DISCUSSION 

To define the orientation of the spin plane in the space, 
it is convenient to introduce the vector normal to the spin 



plane n = [Zi x Z2], where li and I2 are the order param- 
eter components introduced in Eq.([T]). Orientation of the 
spin plane in a magnetic field is determined by the com- 
petition of the anisotropy of magnetic susceptibility and 
the crystal anisotropy. The Zeeman energy and energy of 
magnetic anisotropy for planar spiral spin structure can 
be written as: 

E^-^H^-^^^^{n-Hf + a,nl+a2nl, (2) 

where x\\i X-L arc magnetic susceptibilities for field ap- 
plied parallel and perpendicular to vector n (or, perpen- 
dicular and parallel to the spin plane, correspondingly). 
The constants of anisotropy are positive and ai > 02 > 0; 
X and y axes arc aligned along the a and b axes, corre- 
spondingly. Such choice of parameters ensures at zero 
field planar spiral structure with n \\ c, as it was sug- 
gested in Refl23. 

As it was argued above, the phase transition at Hd 
from the phase I to the phase II is the reorientation tran- 
sition. The observation of spin reorientations for two field 
directions H \\ a, H \\b is possible only for > x±'- the 
spiral structure rotates in strong field to provide n \\ H . 
The critical fields for the H \\ a,b can be obtained by 
minimizing of the energy (Eq. 2) over the direction of 
the vector 

H\\a: 

H.. = .P^, (3) 

V ^11 ^ 

H \\b : 

H.. = .p^, (4) 

V ^11 ~ 

The values of fi^Hd for H \\ a and H \\ b differ 
essentially: 15 T and « 1.8 — 4.2 T respectively (see 
phase diagram on Fig. 7). Thus, we can conclude that 
01/12 « 12 — 70. The anisotropy parameters oi and 
a2 define also the frequencies of electron spin resonance 
(ESR). The spin reorientation transition at ^loHd ~ 16 
T for i? II a was predicted from the ESR data3^. Our 
high field M(II) experiments provide direct experimental 
proofs for this predicted transition at ^J-oHd ~ 15 T. 

The two stage transition from the paramagnetic phase 
to the planar spiral phase is specific for the magnetic sys- 
tems with a strong " easy plane" magnetic anisotropy for 
the vector nM^ This scenario is applicable to the case of 
LiCu202 since ai ^ a2. Thus, LiCu202 can be approx- 
imately considered as a magnetic system with uniaxial 
anisotropy. For these systems at Tc2 < T < Td , only one 
component of the order parameter of the spiral structure 
appears (li || a, see Eq.(IT])) and the other is fluctuating, 
while both components order at T < Tc2- 

Peculiarity of the crystal a-direction is also confirmed 
by the magnetization measurements. The temperature 
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derivatives of the magnetization (Fig|4]) for H \\ b, c 
are similar to each other, and they demonstrate sharp 
anomahes both at Td and at Tc2- In the same time, the 
dM/dT curve for H \\ a changes essentiaUy only at Tc2- 

In the high field range, the second phase transition 
was observed at Hc2- This transition is accompanied by 
the decrease of the susceptibility and is observed for all 
three principal field directions. This fact suggests the 
exchange origin of this transition. At Hc2^ a spiral struc- 
ture in individual chains probably changes to the other 
one. Such a transition from the spiral to collinear spin- 
modulated phase was observed recently in the quasi-one 
dimensional magnet with the same type of frustration 
LiCuV04.— For the quasi-one dimensional systems, this 
unusual magnetic phase was interpreted theoretically as 
a spin-density wave phase^r— . 

For the one dimensional frustrated chain with the ex- 
change constants Ji and J2 of LiCu2 02 (J1/J2 ~ —2) 
corresponding transition is expected, according to Ref|7|, 
at Hc2 ~ The value of -ffsat is not known 

experimentally, but it can be evaluated by extrapolat- 
ing the M{H) dependences to the magnetization value 
equal to gSfiB per Cu^"*" ion. The linear extrapolation of 
M{H) gives the upper estimate for the saturation field 
Hsat ~130 T. Taking into account the 15 percents spin 
reduction in low field range, we can expect, that the sat- 
uration field will be w 15% less than the value derived 
from linear extrapolation: Hsat ~ HO T. This results in 
the evaluation of the Hc2 value of w 20 — 25 T which is 
close to the values observed in experiment. The theoreti- 
cal consideration of the same ID model in Ref. predicts 
a transition to the spin-density wave phase at the value of 
magnetic moment 0.l2Msat, which is also in the reason- 
able agreement with the experimentally observed value 
«0.07M,at (see Fig.©. 

It is unclear how sensitive is the spin-density wave 
phase, predicted theoretically for ID system, to the in- 
terchain exchange interaction J3, which is essential for 
LiCu202 . The effect of the interchain exchange inter- 
action may be responsible for the difference between the 
results of the model calculations and experimental values 



of Hc2- Thus spin-density wave phase can be a plausi- 
ble candidate for magnetic phase realized in LiCu202 at 
H > Hc2- More definite conclusion requires further ex- 
perimental study of this phase. 

Finally, it should be noted that the uniaxial anisotropy 
observed in LiCu202 is unclear at the moment and still 
awaits theoretical consideration. 



VI. CONCLUSIONS 

The low field magnetic phase transition observed on 
the untwinned LiCu202 crystals for the field aligned in 
the (a6)-plane of the crystal can be described as a re- 
orientation transition of the spiral structure. The ob- 
served phase transitions are in agreement with the pres- 
ence of the strong easy (6c)-plane anisotropy for vector 
n of spiral structure in LiCu2 02^. The origin of such 
strong anisotropy in LiCu2 02 at the moment is unknown. 
New high field magnetic transition was observed for the 
field directed along all three principal axes of the crystal. 
The critical field of this magnetic transition is close to 
the value of transition field from a spiral to a spin den- 
sity wave structure predicted theoretically for one dimen- 
sional model with parameters of exchange interactions Ji 
and J2 expected for LiCu202. 
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